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Abstract. The Lindblad equation governs general markovian evolution of the density 
operator in an open quantum system. An expression for the rate of change of the 
Wigner function as a sum of integrals is one of the forms of the Weyl representation 
for this equation. The semiclassical description of the Wigner function in terms of 
chords, each with its classically defined amplitude and phase, is thus inserted in the 
integrals, which leads to an explicit differential equation for the Wigner function. All 
the Lindblad operators are assumed to be represented by smooth phase space functions 
corresponding to classical variables. In the case that these are real, representing 
hermitian operators, the semiclassical Lindblad equation can be integrated. There 
results a simple extension of the unitary evolution of the semiclassical Wigner function, 
which does not affect the phase of each chord contribution, while dampening its 
amplitude. This decreases exponentially, as governed by the time integral of the square 
difference of the Lindblad functions along the classical trajectories of both tips of each 
chord. The decay of the amplitudes is shown to imply diffusion in energy for initial 
states that are nearly pure. Projecting the Wigner function onto an orthogonal position 
or momentum basis, the dampening of long chords emerges as the exponential decay 
of off-diagonal elements of the density matrix. 
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1. Introduction 

Wigner functions, W{x), represent a density operator, p, in the phase space, x = 
{p, q) = (pi, ...,pi, qi, qi), for a system with I degrees of freedom. Given the position 
representation the Weyl- Wigner transformation defines |1[ 

W{x)^ Jdq'{q+y'\j£-^^\q-lq') exp[-'-p-q'y (1.1) 

If we replace within ( p..lD p/(27r?i)' by an arbitrary operator B that acts on the Hilbert 
space of the quantum system , the resulting function, B{x), is then known as the Weyl 
representation, or the Weyl symbol for B. 

The semiclassical approximation for pure states, p = has the familiar WKB 

form 



j 



;i.2) 



where the Sj's are classical actions defined on a classical manifold, i. e. the integral 
of p • dg from some arbitrary position. The corresponding Wigner function, obtained 
by inserting (|1.1| ) in (|T]^) and evaluating the integral by stationary phase, can be 
decomposed into 

W{x)^J2^k{x), (1.3) 
j 

where each term has the form 

Wk{x) = A,{x) cos^. (1.4) 

n 



The sum in ( |1.3|) runs over all chords centred on x that connect appropriate paths 
on the classical manifold manifold that corresponds to For instance, in the case 
where / = 1, this is just the Bohr-quantized energy shell and 5*^ is the area between the 
chord and the shell (plus a semiclassically small Maslov phase) as shown in Figure 1. 
Both semiclassical approximations ( |1.2| ) and ( |1.4| ) break down along caustics, where the 
amplitudes aj, or Ak, display spurious divergences. The caustics of Wigner functions are 
the loci of coalescing chords. It is there nescessary to resort to uniform approximations 
for the integral (1.1), resulting in Wigner functions that have large but finite values 
in the neighbourhood of the caustics. Detailed studies of these are given by Berry 
0, for / = 1, and in for general integrable systems. The various expressions for 
the semiclassical amplitudes in ( |1.4| ), also given by these references, are presented in 
the Appendix. The verification that ( |1.4] ) is consistent with the pure state condition, 
p^ = p, within the semiclassical approximation, irrespective of caustics, was carried out 
in i. 

In order to describe mixed states, 

P = Y.(^n\^n){^n\, (1-5) 
n 

we merely superpose the corresponding semiclassical Wigner functions. This is simpler 
than in the position representation, in which we must multiply the sums for the bras and 



Decoherence of Wigner Functions 



3 




1 



Figure 1. A single cliord, ^, is centred on x if it is close to a convex energy shell in 
the case of one degree of freedom. The phase of the Wigner function is proportional 
to the area S{x) between the chord and the shell, while the ampHtude, A{x), depends 
on both phase space velocities at the tips of ^, as described in the Appendix. 



the kets. Indeed, there are cases where the semiclassical approximation to the mixed 
Wigner function is approached directly, instead of relying on the pure states [Q. Such 
is the case of the spectral Wigner function [^]. This is defined as 

{2nhy ^^^^ ^ ^ _ \ {E-Enf 



W(x;E,e) = 



X) exp 



2e2 



;i.6) 



where VFn(x) is the pure state Wigner function for the n'th eigenstate of a given 
Hamiltonian. If the width of the energy window, e, is classically small, but contains 
many eigenstates, then the semiclassical limit of ( |1.6| ) can also be decomposed into 

(1.7) 



W(x;E,e)^^Wfc(x;E,e), 

k 



where 



Wfc(x;£',e) = A exp 



COS 



1.8) 



Here iSfc(x) is again the area closed off by a chord centred on x, but both its tips are 
now further constrained to lie on the same classical trajectory, within the energy shell 
of energy E. The time taken for the traversal of this trajectory segment is and the 
amplitudes Ak are given in references If an alternative Laurencian energy window 

is used, such as in [§], then the chord contributions decay exponentially, rather than 
with a gaussian factor. It is important to note that (|1.^ ) holds for any kind of classical 
motion, whether integrable, or chaotic to any degree, since it depends only on finite 
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times of the order of fi/e. In the case of chaotic motion, the semiclassical hmit of the 
individual Wn{x) in ( |1.6| ) is not known, but we can still work with the semiclassical 
mixture ( |1.^ ). 

Even initially pure density operators evolve into a mixture if the system is open. 
Thus the Wigner function of a system that attains equilibrium with a heat bath at the 
temperature T is given by the canonical distribution 

W'.W~i:exp(-§)^W. (1.9) 

where Z{T) is the partition function. The semiclassical limit of Wt{x), in the case of 
/ = 1 was obtained by Korsch 0. 

How can we take the external environment into account in the evolution of an 
initially pure state? The influence functional approach to this problem is to construct 
a specific (simplified) model of the environment and its interaction with the system, so 
as to subsequently trace out the external variables [|TU], [TI| . Alternatively, the theory of 
semigroups starts from the mathematical restrictions that the resulting mixed density 
matrix must satisfy. A remarkable simplification follows from the assumption that the 
environment reacts to the system sufficiently fast so as to loose all memory of prior states 
of the system, i.e., the evolution is assumed to be markovian. The density operator must 
then satisfy a linear differential master equation of the following form: 

^ = -{[H,P] + -j:iL,pL] - \L]L,p - \pL]L,). (1.10) 



This theorem represents a culmination of work on open systems [0, |T^, |1^, |T5[ and is 
known as the Lindhlad equation. 

The first term accounts for the unitary evolution of the closed system under its own 
internal Hamiltonian operator, H. The Lindhlad operators, Lj are arbitrary in principle, 
but physically they must be related to the observables which couple the internal system 
to the environment. The simplest example is that of a system that is subjected to an 
external field, so that the interaction depends on its position, but where there is no 
net exchange of energy with the environment. In this case, we can assume that the 
Lj = cjqj, i. e. the Lindblad operators are proportional to the components of the 
position operator. The dependence of the coupling constants on h is already divided out 
1.10|) as discussed in |16[. Once the master equation has been constructed on the 



m 



basis of physical assumptions about the way that the system is affected by the external 
world, the dynamical variables of the environment play no further explicit role in the 
internal evolution. It is important to distinguish the Hamiltonian appearing in ( |1.10| ), 
which drives the internal motion, from the one that defines the pure, or mixed initial 
states in previous equations. If they happen to coincide, then the density operator will 
be constant in the absence of interaction with the environment. 

The purpose of this paper is to analyze the evolution of Wigner functions that 
have initially the semiclassical pure state form ( p..3|) , or the spectral form (|1.7| ), under 
the infiuence of an arbitrary Hamiltonian and arbitrary Lindblad operators. It goes 
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without saying that alternative methods to the decomposition of master equations 
for the density matrix in any orthogonal basis are much more valuable than for the 
Schroedinger equation itself, since the number of states in the truncated basis needs to 
be squared. This is the reason for the popularity of stochastic formalisms reviewed in 
p!6| , but the simple adaptation that needs to be made for the inclusion of hermitian 
Lindblad operators in the semiclassical evolution allows us to work directly with the 
Wigner function and its projections. 

The fact that both the Weyl- Wigner transformation (|1 . 1|) and the Lindblad equation 
( |1.10| ) are hnear with respect to superpositions of density operators allows us to deal 
separately with each chord contribution, Wk{x), or Wk{x;E,e), as if it represented a 
chord density pk- The only added assumption is that both H and the Lj are represented 
by smooth phase space functions in the Weyl representation. This is the usual case 
for an isolated system, where the observable H is given in phase space by a function 
that is at least very close to its classical limit. Though lindbladian dynamics are less 
familiar, the more commonly considered examples 0] couple the internal system to 
the environment through observables, represented by real phase space functions, or, for 
instance, the annihilation operator of the harmonic oscillator, which is represented by 
a complex (but smooth) function. 

No attempt will be made to analyse the accuracy of the approximate Markovian 
evolution in a concrete physical situation. This is merely an extension of the usual 
procedure followed in ordinary hamiltonian quantum mechanics: Though no system 
is completely isolated, it is worthwhile to build general quantum theories under the 
assumption of isolation, with the eventual option of later including external interactions 
||. Likewise, we can now analyse the semiclassical limit of lindbladian dynamics, 
which in some cases should be a preliminary step in the study of more complex non- 
markovian interactions with the environment. Evidently, general lindbladian dynamics 
(semigroups) also encompass the ordinary (unitary group) evolution of closed systems. 
It will be shown that the effect of the environment is easily accommodated within the 
semiclassical evolution in the case of hermitian Lindblad operators. The analysis of 
the evolution of a measure of decoherence, the linear entropy is part of the subject 



of a companion paper ||T8[, though the numerical examples there use finite time step 
equivalents of lindbladian evolution. 

It is important to establish distinctions and connections between this paper and 
previous work on the semiclassical limit of open systems. Perhaps this is the first time 
that the theory is developed entirely in phase space. Although Caldeira and Leggett 
|[TT| already make some use of the Wigner function, it appears only in an auxiliary 
role. This is also the case of Grossmann |T9[ who discusses the semiclassical limit of the 
influence functional approach. A remarkable paper by Strunz comes much closer 
to the present work. A path integral is there developed, but, unlike usual influence 
functionals, it is constructed for markovian Lindblad evolution. The propagator is 

§ For a recent discussion of non- markovian stochastic treatment, see e. g p7|. 
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defined in the position representation, using the Wigner function. Its semiclassical limit 
is derived, but this does not describe immediately the propagation of the semiclassical 
Wigner functions themselves. Indeed, this is a tricky problem even for pure Hamiltonian 
evolution, that has been fully understood only recently |]21], The important point 



is that Strunz defines a decoherence influence functional for pairs of arbitrary paths in 
phase space that is here shown to account for the dampening of the amplitude of each 
chord of the semiclassical Wigner function. The same decaying factor is easily seen 
to affect the off-diagonal elements of the semiclassical density matrix in the position 
representation. 

Thus, the next section reviews the unitary evolution of semiclassical Wigner 
functions and discusses the purity of the states. This is followed in section 3 by the 
derivation of the Lindblad equation for the Wigner function without the hamiltonian 
term. The equation is integrated in the case of hermitian Lindblad operators. Section 
4 combines the hamiltonian and lindbladian evolutions of both preceding sections to 
produce the semiclassical description of decoherence of an evolving Wigner function 
over a finite time interval. Energy diffusion of a stationary Wigner function in an 
open system is discussed in section 5. Section 6 shows how the inverse Weyl- Wigner 
transformation leads to a semiclassical approximation in the position representation 
that is again described by the decoherence distance functional. The delicate issue of 
normalization is treated in the Appendix. 

2. Review of unitary evolution 



The Moyal Formula 

^W{x) ^ {H{x), Wix)} + 0(fi2), (2.1) 

where {,} denotes the classical Poisson Bracket, suggests that the semiclassical evolution 
follows the classical trajectory for the argument of the Wigner function. This is indeed 
the exact result for quadratic Hamiltonians and also holds approximately if W{x) is a 
smooth function within the scale of H. For pure semiclassical states such as ( |1.4|) , or 
for the mixture ( |1.8|) this is certainly not the case. However, it has been shown that 
the evolution then depends only on the orbits of both tips of each chord centred on x 



21| , |22| . Since neither of the two previous treatments of this problem relied directly on 
the differential equation for pt, i.e., the unitary part of (|1.1CI|) , it is a useful exercise to 
rederive the hamiltonian evolution in this way, preliminary to the inclusion of the new 
terms in the master equation. 

The main ingredient is the integral form for the product of two operators B2B1 in 



the Weyl representation, due to Berezin ||2^ 

21 



ft 



(2.2) 



B2Bi{x) = (^-^^ J dx2dxi i?2(a;2)-Bi(xi) exp 
which was subsequently extended to higher order products [^, 0] . Here, 

A(x, Xi, X2) = 2(x A Xi + Xi A X2 + X2 A x), (2.3) 
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where the skew product a Ab = (Ja) ■ b, with 














[2A) 



[2.5) 



the matrix that exchanges positions and momenta in Hamilton's equations: 
dx jdH 
dt dx ' 

Each term in (|2.3|) is the (symplectic) area of the parallelogram defined by the 
corresponding vectors (i.e. just the area itself for / = 1). Much more useful is to consider 
A as the symplectic area of the unique triangle that has {x,Xi,X2) as its midpoints. 
This is the circumscribed triangle, rather than the usual inscribed triangle defined by its 
corners. 

The equation that governs unitary evolution of the Wigner function is thus 



d 2 / 1 /■ 

—W(x) = - — / dx2dxi H(x2)W(xi) sin 
at h Vvrfiy J 



— {x,Xi,X2} 
n 



[2.6] 



The evolution of semiclassical Wigner functions can be decomposed into separate chord 
contributions. Hence we insert ( |1.4D , or ( |l.(j| ) into ( |2.(j| ) and evaluate the integral within 
the stationary phase approximation [[]]. Expressing the sines and cosines as the sum of 
exponentials, there will be four integrals for each chord of the Wigner function: 



d_ 

di 



W,{x] 



1 

'2h 



1 



21 



E 

± 



da:;ida;2 H{x2)Ak{xi) 



exp 



-[A±S,(xi)] 
ft 



exp 



-\/\±Su{x,)] 



(2.7) 



The stationary phase points are now obtained from the rule that the derivative, d/dx, 
of a symplectic area, bounded by a chord ^, is just 0. Hence, from 

{d/dxi)[^±Sk{xi)\ = -2 J (x2-x)± Ja(xi) = 

{d/dx2)[A ± Sk{x2)] = 2 J (xi - x) = 0, (2.8) 

it follows that the circumscribed triangle. A, collapses at the stationary configuration 
for which xi = x and C,k{xi) = ±2(x2 — x), i. e. the Hamiltonian must be evaluated at 
the stationary points, which coincide with the tips of the chord centred on x. At these 
points, the hessian determinant of the phase is just 

2J 



det 



dxidx2 

2J 







(2.9) 



and the signature of the hessian matrix (the number of positive, minus the number of 
negative eigenvalues) is zero. To obtain this result, note that it obviously holds for 
S = and we can always connect the eigenvalues of (|2^ ) continuously to this particular 
case (avoiding caustics). Since the eigenvalues remain finite along any path and their 
product is unity, the signature remains zero. 



II The following development will be made for the pure Wigner function (1.4). This holds also for the 
spectral Wigner function (|l.S|), by including the exponential factor in the amplitude. 
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Combining all these ingredients, we obtain the stationary phase evaluation of (127 



as 




H\x-^ 



Adx) sin 



Skix) 



(2.10) 



which can be immediately compared with the straightforward differentiation of (|T7 
The dominant term is 



di 



sm 



Sk{x) 



which agrees with ( p.lO|) , provided that 



dt 



H{x + 



H 



X 



(2.11) 

(2.12) 
is just the 



But this form of the Hamilton- Jacobi equation, derived explicitly in p2 
derivative of the action in It differs from Marinov's Hamilton- Jacobi equation 

which depends on just a single tip of the chord 



Thus, we find that the stationary phase evaluation captures the dominant term 
of the unitary evolution of the Wigner function, while neglecting the variation in the 
amplitude Ak{x). This propagation depends on the classical flow, through the orbits 
of both tips of each chord, ^k{x). (Only if H{x) is quadratic, will the difference of 
the hamiltonian vector fields at the chord tips in (|2.1CI|) coincide with JdH/dx at the 
centre.) Of course, we can also evolve the entire classical manifold and then reconstruct 
the Wigner function at the centre of the evolving chords as described by Berry and Balazs 
p7| , resulting in a relatively slow change of amplitude ^ that is not captured by 

To conclude this section, we note that both pure and mixed states can undergo 
unitary evolution, but that, in the former case, the purity of an initial state is not 
destroyed, i. e. we must have fP' = p for all time if this is true initially. In the Weyl 
representation, the product rule ( ^^ leads to 

2 \^ f 

dxidx2 W{xi)W{x2) exp 



W{x) 



n 



A{x,Xi,X2) 



(2.13) 



which was shown to hold for general integrable systems within the stationary phase 
approximation [Q. It is quite remarkable that this result is entirely impervious to the 
multiplicity of chords or their coalescense at caustics, which are washed out by the 
integration. 

An important consequence of the fact that pure state Wigner functions satisfy the 
pure state condition is that we can take the trace of both sides to obtain 

tr p = JdxW{x) = (27rfi)' J dx [W{x)f. (2.14) 

It happens that the integral involving the square of the Wigner function does indeed 
equal unity, within the stationary phase approximation, while the first integral does not, 
though its semiclassical value is independent of Planck's constant. These results were 
first obtained by Berry [@], for / = 1, and are extended to general integrable systems 
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in the Appendix. It may be said that the pure state Wigner functions are indirectly 
normalized. 



3. Pure lindbladian evolution 



The Weyl representation of the remaining terms of the master equation ( |1.1CI| ) requires 
the integral form for the product of three operators B-^B2Bi. Sometimes, as in it 
is necessary to use the full generalization of ( |2.2|) , integrating over three phase space 
variables with a phase that is a circumscribed quadrilateral . However, in this case it 
is possible to use the following shortcut also provided by 0: 



21 



B:iB2Bi{x) = { — \ [ dx2dxi Bs{x2- xi + x)B2{x2)Bi{xi)exp -A{x,Xi,X2) ■ (3.1) 

Thus, the lindbladian evolution for a single chord of the semiclassical Wigner function 
( |1.4D becomes 

(9 ~ 1 / 1 \ 

^Wkix)^(—) V 



2h Vvrfi 



|y dx2dxi Lj{x2 — Xi + x)Ak{x2)Lj{xi)* exp 
-^jdx2dx, L,{X2 -x, + xrL,{x2)A,{x,) exp 



{A{x,Xi,X2) ± Sk{x2)) 



- {A{x,Xi,X2) ± Skixi)) 
n 



(3.2) 



if we neglect the effect of the Hamiltonian. Here, 3? stands for the real part and we note 
that an adjoint operator W is represented by the complex conjugate function B{x)* . 
Presuming that the Lindblad operators are represented by smooth phase space functions, 
the stationary phase evaluation of ( |3.2| ) proceeds essentially as in the previous section. 
Collecting the terms, we obtain 



Ljix + y)i^i(a; - y)* exp (^^Skix) 



\L,{x+^p' + \L,{x-^p^ 



cos 



Sk{x) 



(3.3) 



From here on the discussion will focus on the case of hermitian Lindblad operators, 
i.e. it will be assumed that the system is only coupled to the environment by observables. 
Then the master equation for Wk{x) simplifies to 

2 



^W,{x) ^ --Akix)cos^^ V 

dt ^ ' 2n ^ ' n ^ 



(3.4) 



In contrast to the case of Hamiltonian evolution, we see that hermitian Lindblad 
operators do not affect the phase of the semiclassical Wigner function, but only its 
smooth amplitude function. Indeed, if we assume that the chords, ^k{x), and the 
actions, Sk{x), remain constant at each point, x, we can integrate the equations for 
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the amplitudes, 



•9 , / N 



3 



to obtain 



Ak{x,t) = exp 



■-E 



Lj{x 



2 ' 



2 ' 



0). 



(3.5) 



(3.6) 



Thus, it is easily verified that the semiclassical Wigner functions ( |1.4| ), now taken with 
the time dependent amplitude (|3.(i|), i. e. 

Wk{x, t) = Akix, t) cos [Skix)/H], (3.7) 

satisfy the evolution equation (|3.4|) and, of course, an equivalent result holds for each 
chord of the spectral Wigner function, Wkix; E, e). 

Evidently the extreme case of pure lindbladian evolution will be relevant only if 
the coupling constants implicit in the Lj are large enough for the internal motion of the 
system to be effectively frozen within the decay time for the amplitudes. The Lindblad 
operator with the largest difference between the tips of a chord dominates the decay and 
the corresponding time is inversely proportional to the square of the coupling constant. 

The simplest case is that of a single operator L. Clearly the Wigner function 
remains constant if the state is initially an eigenstate of L, since the tips of all chords 
then lie on the same manifold L{x) = constant. This result holds for all the cases 
discussed in the Introduction: For I = 1, this manifold is just a phase space curve. For a 
general integrable state, the chord tips lie on an /-dimensional torus, within the (2/ — 1)- 
dimensional surface, on which L{x) = constant 0. The chords of a spectral Wigner 
function ( |1.8| ) also lie on the single surface H{x) =constant. Hence, if H{x) = L{x), the 
Wigner function is also invariant. 

Generally, for H{x) ^ L{x), it must be that L{x+) ^ L{x-) for most chords, 
^ = a;+ — x_ . Evidently, short chords are more stable against decay due to arbitrary 
Lindblad operators, but, for specific operators, there exist long chords that will also 
be selected for survival. For instance, suppose that / = 1 and L{x) = cq, the position. 
Consider a given centre, x, with a vertical chord, C,{x), as in Figure 2, then the amplitude 
does not decay. In this case, the full Wigner function condenses onto the classical curve 
and onto the line of midpoints of the vertical chords that bisects the classical curve. 
However, we shall see in the following section that this strange effect only arises in 
the complete absence of the hamiltonian term in the evolution equation: No matter 
how strong the coupling to the environment, the Hamiltonian will eventually rotate the 
chords that are then erased by the Lindblad operator. For a general position coupling 
to the environment, with / > 1, i.e. L = (cigi, ...,ciqi), the decay rate will depend on 
the distance: [ci(gi+ - gi_)^ + ... + Ci{qi+ - qi^f]^^'^. 

Evidently, the above picture is entirely consistent with previous results on 
decoherence [111], |16|, EH]. The general initial pure state may be considered as a 



superposition of the eigenstates of the Lindblad operator. The effect of "continuous 




Figure 2. Chords of the energy shell with the same q± can be short, ^i, or long, ^2- 
Both will be equally dampened by the Lindblad operator L = cq. After a long time, 
the Wigner function condenses onto the energy shell and the line of midpoints of the 
vertical chords. 



measurement" of this observable is to erode nondiagonal elements of the density matrix 
in this representation If we merely place a pair of coherent states at both tips of 
a given chord, we obtain a Wigner function that is a gaussian at x± = x ± with a 
width of 0{\/h). The interference between these two states is represented by another 
gaussian centred on x itself but this is modulated by narrow fringes, of wavelength 
0(/i). This central gaussian decays exponentially just as the general semiclassical Wigner 
function that we have been studying when a Lindblad operator is "turned on", unless 
L{x+) = L(x_). 

An initial semiclassical Wigner function of the type ( p. .41) , with the amplitudes 
specified by (|A.7|) , becomes mixed due the lindbladian decay of the amplitude (|3.6|) . 
The way that this decoherence is connected with energy diffusion, within the framework 
of the spectral Wigner function, will be discussed in Section 5. An important issue is the 
normalization of the mixed Wigner function that results from lindbladian evolution. It 
is shown in the Appendix that the pure state itself is only indirectly normalized through 
tr/5^, but this is now a measure of the decoherence of the mixed state and should not 
equal tr p, once the environment is switched on. However, we find that the semiclassical 
evaluation of tr p, though it is not unity, is stable with respect to the lindbladian decay 
of amplitude, as well as with fi, taken as a small parameter. Therefore, we can safely 
attribute the semiclassical change in tr as entirely due to decoherence. 
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4. General evolution 



We are now in a position to combine the unitary evolution of a closed quantum system 
with the decoherence that results from its coupling to the environment, as modeled by 
Lindblad operators. Again we shall treat each chord contribution independently, just 
as in the two previous sections. The differential equation corresponding to the master 
equation ( |1.10D is just the combination of (|2.1CI| ) with ( p. 41) in the case of hermitian 
Lindblad operators, but it is hard to integrate this directly. Instead, the evolution can 
be taken as the result of alternate sequences of unitary and pure lindbladian evolutions in 
the limit of arbitrarily small time steps. This approach has been used for the coupling 
of two different unitary maps by Berry et al. pOf. The recent study of a unitary 
map alternating with a finite step lindbladian evolution by Bianucci et al. pOl is even 



more to the point. The important issue is that the exponential damping factor in the 
amplitude remains constant during the unitary evolution described in section 2 through 
the classical motion of each chord and its centre. Alternatively, the chord is then frozen 
during the interval in which the system is governed only by the Lindblad operators. 

Let us define Vk{x,t) as the semiclassical contribution of each chord to the Wigner 
function in the absence of the Lindblad operators, i. e. Vk{x,t) describes purely unitary 
evolution. Splitting the period of evolution, t, into N half-steps of hamiltonian evolution 
and half-steps of lindbladian evolution, we obtain 



N 



Wk{x,t) = Vkix,t) Yl exp 



n=l 



t 



2nN 



E 



Lj Xn + 



^k{Xr 



Xr. 



ik{Xn) 



(4.1) 



Naturally, it has been necessary to multiply the Hamiltonian by a factor of 2 in this 
formula, since it only acts during half of each interval, and the coupling constants in the 
Lindblad operators are likewise increased by a factor of \/2. The final dampening of the 
amplitude at x = x^v depends on both the history of the centres, and that of their 
chords, ^k{,Xn)- In other words, the total decay depends on the classical trajectories 
that arrive at Xk,N± = a;Ar ± ikix^)- Each factor in the product is then calculated at 
the centres, Xk,n = {xk,n+ + Xfc,„_)/2, with the chords, C.k,n = Xk,n+ - Xk,n-, during the 
half-intervals in which the hamiltonian is switched off. (Note that the sequence of Xk^nS 
depends on the final chord ^k, as well as on x.) 

Taking the limit A^ ^ oo, we obtain the continuous evolution. 



Wkix,t) = Vk{x,t) exp 



1 

2fi 



/ dt'\Ljixk+it'))-Ljixk-ii') 



.(4.2) 



Therefore the full evolution of the Wigner function in an open system combines the 
unitary evolution for each chord, obtained by the classical motion of the chord tips, 
Xk+(t) and Xk-it), with the exponential decay of the amplitude, characteristc of pure 
lindbladian decoherence. But now this dampening depends on the chord history. If 
there is more than one chord centred on x, each contribution to the evolution of the 
Wigner function depends on the particular motion of both tips of each respective chord. 
Generally the classical motion of the centre of a chord does not coincide with the 
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midpoint of the motions of the tips, unless the motion is hnear (i. e. the internal 
Hamiltonian is quadratic). Even though the development has focussed on the Wigner 
function of integrable systems, exactly the same form of general lindbladian evolution 
also holds for each chord contribution, Wfc, to the spectral Wigner function ( |1.8| ). The 



normalization of ( ^.2| ) is discussed in the Appendix. 

It is now feasible to make a partial connection to the path integral developed 
by Strunz. Indeed, the phase space influence functional in |^ has exactly the same 
form as the exponential in ( [4.2|) in the case of hermitian Lindblad operators, but 
this is there taken over all possible pairs of paths, which are not restricted to being 
classical trajectories. Even so, working within quite diverse frameworks, both treatments 
conclude that the effect of hermitian Lindblad operators is encapsuled by the decoherence 
distance functional: 



Dt[x+{t'),x4t')] 



1 



Y: f dt' \L,{x^{t') - L,{x.{t'))\ 

■ JO 



(4.3) 



Though Strunz also discusses the semiclassical limit, this is treated within the context of 
an effective Hamiltonian, H{p, 5') — ^ l-^jP? leading in general to complex trajectories. 
However, in the specific example of the harmonic oscillator with the Lindblad operator 
L = cq, the classical paths in that contribute to the propagator in the position 



representation turn out to be real. 

The motion of the chords is important even for the Wigner function corresponding 
to an eigenstate of the internal Hamiltonian, i. e. even if Vk{x, t) is stationary. The tips 
of each chord then evolve simultaneously along trajectories constrained to the invariant 
manifold (just the energy curve if Z = 1) that defines the semiclassical Wigner function. 
Thus, the dampening of the amplitude is only small for short chords, so that the line of 
midpoints of vertical chords in Figure 2 is erased if the internal Hamiltonian is taken into 
account. The same is true for spectral Wigner functions of the internal Hamiltonian. 
Thus, unless all the {Lj, H} = 0, only the classical region of the Wigner function close to 
the energy shell survives the decoherence due to the external environment. This general 
result is compatible with a previous specific deduction for the harmonic oscillator 



5. Energy diffusion 

Let us now consider the ultimate effect of coupling a closed system to the environment 
through Lindblad operators. The focus will be on an initially mixed state within a 
narrow energy window, corresponding to a spectral Wigner function of the internal 
Hamiltonian. Of course, this also includes a pure state, in the limit where the window 
width tends to zero, for the simplest case where / = 1. The amplitudes of the various 
chord contributions to W(x; e) in (|1.7|) decay exponentially with {Dt[x+{t'),x_(t')]y 
and the trajectories x±{t') are constrained to the stationary energy shell. Hence, only 
short chords with centres very close to the shell will give appreciable contributions to 
this Wigner function after a finite interval of time. 
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There are only two possibilities concerning the trajectories that connect the tips 
of a very small chord, ^k'- Either this is the unique short trajectory that takes a small 
time, Tfc, between the two tips, or else the trajectory winds near a periodic orbit ||^, |^. 
For / = 1, the energy shell is itself a single periodic orbit, so that the difference between 
both types of contributions concerns only the number of repetitions around the periodic 
orbit (which is zero for the short orbit). Both types of contribution will arise for all 
centres close to the shell in this simple case. For / > 1, the long orbits that contribute 
to the Wigner function with a small ^fc(x) will be associated to neighbouring periodic 
orbits with varying actions, as x is moved along the {21 — l)-dimensional shell. Thus, 
for any number of degrees of freedom, we can separate 

W(x,t;E,e) — >yVr{x,t;E,e) + moXx,t;E,e), (5.1) 

t — *oo 

where we lump all the periodic orbit contributions into the second term. It is important 
to distinguish here the duration, t, of the evolution of the open system, from the time, 
r, that the classical trajectory takes to travel between the tips of the short chord. Thus, 
r is a function of x, for a fixed energy surface and is completely independent of t. 
It is easy to estimate the short orbit contribution, because 

^T-(x) ^ XT = tJ dH/dx, (5.2) 

to first order in r . Hence, we can approximate the difference of the Lindblad functions 
in terms of a Poisson bracket: 

|L,(x + e./2) - L,ix - e./2)p ^ t'\{H, L,}p. (5.3) 

Let us now assume that the closed classical system is ergodic. Then the time 
average of \{H, Lj}\'^ for a trajectory on the energy shell exists and generally equals the 
average over the shell. But, for a short chord, the classical trajectory of the midpoint 
x{t) ^ [x+(t) + x_(t)]/2, so that, considering now r to be an average time for the 
traversal from x_(t) to x^{t), we can approximate the decoherence distance functional 
(14. 31) so as to obtain, 



Wrix, t; E, e) ^ V^(x; E, e) exp 



(5.4) 



where we recall that in this case Vr, the Wigner function for the closed system, is 
stationary. Thus the short orbit contribution to the spectral Wigner function, after a 
time t of contact with the environment, has the same form as the initial spectral Wigner 
function (|1.8| ). However, the width of the energy window grows as 



ht 



m? = ^'+i^i: \{H,L,}\\ (5.5) 

i.e. the internal energy of the system spreads diffusively. This will also be the case for 
an initially pure state if Z = 1. The absence of any net drift in the semiclassical motion 
for a mixed state under the action of hermitian Lindblad operators is consistent with 



previous work on the Lindblad master equation . 
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If / > 1 a pure initial state will surely also mix with neighbouring eigenstates of 
the other variables that together define its invariant torus, but this effect is harder to 
retrieve. If the initial spectral Wigner function is defined for a different observable, 7i, 
than the one that drives the internal motion of the system, the surface 7i(x) =constant, 
which is the locus of the tips of all chords, will no longer be stationary. In this case, we 
can no longer invoke ergodicity to estimate the growth of the Tij-energy window, though 
a similar argument can be used if the movement of the energy shell is adiabatically slow 
as compared with the time that the orbits inside the shell require to cover it. Of course, 
this approximately adiabatic scenario is more likely to arise for / = 1. In any case the 
decay of the amplitude of the very short chords that survive the decoherence after a 
finite time will still be ascribable to a net growth of the energy window. 

If Wt were the only term in the decomposition ( ^.1|) of the spectral Wigner function, 
the effect of the Lindblad operators in the semiclassical evolution would be reduced to 
merely mixing the state in a widening energy window. That is, the density matrix 
would be diagonal in the eigenbasis of the operator used in the definition of the initial 
state. The periodic orbit terms in ( |5.1| ) modify this picture in a way that requires 
further analysis. In other words, scars of periodic orbits will play a significant role in 
the semiclassical theory of decoherence, since their contribution to the evolved Wigner 
function is not merely the same as in a spectral Wigner function that has a very narrow 
energy window: Rather than being dampened by their period, the scars only decay as a 
funcion of the length of the chord, as measured by the decoherence distance functional. 



6. Position representation 



Instead of the basis of eigenstates of the internal Hamiltonian, let us consider another 
orthogonal basis, such as the position representation. Then even the density matrix 
for an initially pure state is viewed as a (coherent) superposition of a large number of 
basis states. In the case of Lindblad operators that were just functions of the positions, 
q, it could have been easier to obtain the evolution of the semiclassical density matrix 
directly from the the position representation of the Lindblad equation (|1.10|) . For general 
hermitian functions of p and q it is better to proceed from the inverse Fourier transform 
to (O): 

q+ + q- 



{q+\p\Q-) = J dp W (^p, ^^"Y^) 



exp 



:P ■ {q^ 



(6.1) 



Sufficiently far from diagonality, i.e. for q+ ^ q-, we need only be concerned with 
the stationary phase evaluation in the oscillatory interior of the Wigner function, so that 
the stationary phase condition is fulfilled where the chord centred at x = (g+ + g_)/2 
has the component ^q{x) = q+ — q-. The lindbladian dampening of the stationary chord 
can be factored out of the integral, so that the density matrix of the evolved system is 



expressed in terms of the decoherence distance functional ( |4.3| ) as 



{q+\p\q-) ^ 



isj+{q+) - Sj_(g_)) 
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exp {-^{A[(g+,P,+)(t'), iQ-.Pj-m]}'} ■ (6.2) 

Here the lable j± specifies the different momenta, pj±, for the torus at the positions 
and the actions Sj± are the usual actions in ( |1.2| ), measuring areas from the g-axes. 

In the hmit as g+ g+, we also obtain a short chord contribution from the 
classical region where the semiclassical approximation breaks down. However, the 
extrapolation of ( |6.2|) for short chords gives the same result as follows from the simpler 
approximation for the Wigner function 0, W{x) = (27r)^'5(/(x) —X) for the eigenstate 
of an integrable Hamiltonian corresponding to the action variables I{x) = I. Indeed, 
the correlation function given in follows from the short chord limit of (|6.2| ) by 
taking (sj+ — Sj_) — > pj ■ (g+ — g_) as g+ g_. Of course, the decoherence distance 
functional vanishes in this limit, so that this contribution becomes insensitive to the 
Lindblad operators. We should note that there will generally also be long vertical 
chords, such as shown in Figure 2, contributing to the density matrix. The projection 
of the semiclassical Wigner function onto position space is a good example of the nice 
behaviour of this approximation within integrals, despite the presence of caustics. 

The contribution of long chords to the density matrix of a mixture of energy 
eigenstates in a classically narrow energy window, corresponding to the spectral Wigner 
function, has exactly the same form as ( |6.2| ). The decoherence distance functional 
dampens the contributions from the various trajectories that satisfy the variational 
principle for fixed energy between q_ and g+, so that, if we lable these by their initial 
and final momenta, pj±, the stationary action is just Sj+ — sj-. The geometry of the 
orbits and the semiclassical amplitudes of their contributions are thoroughly discussed 
by Toscano and Lewenkopf in a recent paper They also provide the limit for the 
contribution of short chords as g+ ^ q-, which tends to Berry's result for / > 1 

oc (6-3) 

[p{q)\q+-q-\/f^Y^^-^ 

where J^, is a Bessel function. This result follows from the approximation of the Wigner 
function as a delta function over the energy shell. Of course, there is no decoherence 
dampening for the null-chord contribution. 

The semiclassical contributions from long chords to are obtained by 

evaluating a similar integral to ( |0| ) (see e.g. 0]) by stationary phase, with analogous 
results to ( [6.21 ). However, the limiting form for the contributions of small chords to the 
projection of the spectral Wigner function does not generally have the same form as 
(p. 3D, because this was obtained under the assumption of a simple p"^ dependence of the 
Hamiltonian. 

The same general picture holds for any basis resulting from a linear canonical 
transformation, x x', of the classical phase space, because of the well known 
symplectic invariance of the Weyl- Wigner formalism (see e.g. [0). The density matrix in 
the g'-representation follows from the transformation (|6.1| ), after performing the classical 
change of variables. For each of these representations the decay of the off-diagonal 
elements of the density matrix {q\\p\q'_) is determined by the decoherence distance 
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functional ( [4.3| ) evaluated at the tips of the chords that project down precisely onto 
and q'_. This delicate quantum effect quickly drives the density towards a diagonal 
mixture of the same initial set of basis states which were superposed coherently to form 
the initial state, if this was an eigenstate of the internal Hamiltonian, because of the 
absence of drift for hermitian Lindblad operators. The ultimate effect of the classical- 
like energy diffusion considered in the previous section will be manifest only on a larger 
time scale. Again, it should be pointed out that this qualitative scenario is well konwn 



see e.g. [|ri|, |16[ ^), the novelty being the complete model-independence in which 



these general semiclassical results are here obtained. 
7. Conclusion 

The evolution of quantum systems that are weakly affected by the external environment 
can be studied within the framework of the Lindblad master equation If the 



Lindblad operators are all hermitian, the classical motion that supports semiclassical 
evolution is unaltered by an external environment, the only effect being to dampen 
the amplitude of each chord contribution to the semiclassical Wigner function. In the 
present context, this result embodies the known fact that hermitian Lindblad operators 



produce no net drift [|1^]. The decay of the semiclassical amplitudes was interpreted in 
section 5 as a diffusion of the energy of the initial (nearly pure) state, which is a more 
obvious effect of opening the system than the delicate loss of quantum coherence that 
immediately manifests itself as a loss of amplitude in the highly oscilatory region of 
the Wigner function. An intuitive explanation for this decay was recently put forward 



by Zurek in the restricted framework of interference fringes in the Wigner function 
for a superposition of coherent states. These narrow oscillations arise near the centre 
of the chord separating a pair of gaussians in phase space that would represent each 
single coherent state. Averaging over the juggling of the coherent states, effected 
by the external environment, destroys the interfernce pattern. This same intuitive 
interpretation can now be generalized to the slight motion of the whole classical manifold 
that supports extended semiclassical quantum states. 

General Lindblad operators also account for energy loss. This is in principle 
contained in the differential equation ( |3.3| ) for the Wigner function, though no immediate 
separation of phase and amplitude is then possible. The path integral of Strunz 



does lead to a semiclassical approximation in terms of complex trajectories. However, 
this is derived in the position representation and even if its Wigner- Weyl transformation 
( p..l[ ) were obtained, it would still not be a clear route to the semiclassical evolution of 
Wigner functions. The reason is that, even for unitary evolution, the propagation of 
the Wigner function depends on the initial chord, through the orbit of both its tips. 
Thus, the semiclassical propagator of Wigner functions is paradoxically irrelevant to 
their motion, as discussed in |^|. In short, the general semiclassical description of the 



evolution of Wigner functions is undoubtedly more involved than the present restriction 
to hermitian Lindblad operators, i. e. observables. 
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It is only in this case that the decoherence distance functional ( [4.3| ), defined by 
Strunz for arbitrary phase space paths, localizes on the real classical trajectories with 
end points at the tips of each chord of the Wigner function. The simplicity of the 
Weyl representation of observables as real smooth phase space functions then leads 
by direct integration to simple expressions for the density matrix in the position, or 
the momentum representation, in terms of the same decay factor that accounts for 
the decoherence of Wigner functions. The decay of large chords then appears as the 
dampening of off-diagonal elements in any representation. 

Of course, one should be wary of the extrapolations that are required to synthesize 
such a simple picture. The description of a spectral Wigner function corresponding 
to a mixture of states in an energy window is only supported semiclassically by a 
single energy shell if the width of the window is classically small. Eventually, energy 
diffusion will widen the window beyond the limit of validity of this assumption. But, 
anyhow, the basic hypothesis in deriving the whole semiclassical theory, that only 
the oscillatory complex exponential terms in the integrals (|2.7|) and (|3.2| ) were rapidly 
varying in phase space, eventually breaks down. Indeed, for large evolution times the 
amplitude functions were found to become sharply peaked along the classical manifold 
that supports the semiclassical state. Finally, one can add the further criticism that the 
simple semiclassical approximation, based on independent chords, breaks down in the 
neighbourhood of the classical manifold on to which the Lindblad operators concentrate 
the Wigner function. There, the pure states have been shown to be described by uniform 
approximations 0, Q and ultimately their decoherent evolution should be ascertained. 

Even so, the simplicity of a broad unified view over a wide range of diverse 
phenomena is a rare enough asset to be thrown away easily! For instance, one is familiar 
with the many uses of the Gutzwiller trace formula within various branches of 
" semiclassics" , in spite of the many conceptual problems with its implementation. In the 
case of quantum spectra, it is also possible to find improved approximations that smooth 
over spurious singularities ( see e.g. fS^), but at a cost. The indirect normalization of 
the square of the pure state Wigner function carried out in the Appendix illustrates 
the uncanny selfconsistency of the simple semiclassical approximation, which previously 
passed the more stringent test of the full pure state condition 0. The difficulties with 
chord multiplicities and caustics simply vanish within such integrals. The companion 
paper [|18[ shows, among other things, how the simple semiclassical Wigner function here 
developed can be used to calculate the growth of linear entropy, an important measure 
of decoherence. Indeed, the stability of the direct normalization integral (also in the 
Appendix) throughout the action of the Lindblad operators, allows us to single out the 
decay of the Wigner function amplitudes as a pure decoherence effect resulting from the 
contact of the internal system with the external environment. 
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Appendix A. Normalization 

The normalization of pure state semiclassical Wigner functions is a delicate matter. 
Even though Berry treated the case of one degree of freedom, it is worthwhile to 
rederive the integrals in a manner that is also valid for tori of higher dimensions. The 
main idea, already employed in |^ , is to transform the integration variables from centre 
coordinates, x, to the canonical angles (0_, 9j^) at the tip of each chord on the quantized 
torus that is specified by its constant action variables, I± = T. 

Let the points near each chord tip, x±, take on local coordinates X±, i. e., phase 
space points are given locally by x = x± + X±. Then one can also define local action 



Decoherence of Wigner Functions 



21 



and angle variables as I±{X±) = I{x) and 9±{X±] 
5X± 



d9± ^ dl± 
Hence, for chord tips constrained to I± 



6{x), so that near each chord tip, 

(A.l) 



X, the Jacobian for the change of variables. 



, —^6x= (X+ + X_)/2, is obtained from 



6x 



36, 



59+ + 



dx^ 



66. 



so that. 



det 



d{p, q \ 



det 



1 dp+ 

2 de+ 

ldq+ 

2 ae+ 



1 dp- 

2 86- 

1 dq- 

2 (90_ 



(A.2) 



(A.3) 



Each element in (|A.3| ) represents an (/ x /)-dimensional block matrix and each of its 
elements may be interpreted as a component, pi, ...,pi, or gi, qi, of the velocity vector 
x± for the various choices of actions Ji(x), Ii{x) as the Hamiltonian. In the case that 
/ = 1, these are proportinal to the hamiltonian velocity vectors shown in Figure 1. Thus 
the Jacobian can be rewritten in the form 



det 



d{p, q 



d{e. 



,0+) 



i) Ideti+Ai_| = (l 

which is obvious in the case that / = 1. In the general case this is obtained by the same 
manipulations used in to prove the complementary relation: 
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det{/_,/+}|. 



(A.4) 



det{/_,/+}| 



det 



(A.5) 



It is important to note that the matrix of Poisson brackets, {/_,/+}, is only / x /, 
whereas the Jacobians in both ( [A .41 ) and ( |A.5| ) have (2/) x (2/) dimensions. 

There is just one chord for each pair (^_, 6*+), but the same chord results from the 
interchange of 6- with 6+ , so that the overall transformation of the integration variables 
is 



y ■ dx = 2-(^'+^) ■ I det{I_, I+}| d9_d9^ 



(A.6) 



In passing from integrals over centres to integrals over pairs of angles, all problems 
with chord multiplicities and even caustics disappear. The chord densities really take 
on an independent existence, so that from this new point of view, they may merely 
happen to share the same centre. Whereas in the original semiclassical Wigner function 
a pair of chords with the same centre coalesce as the centre moves on to a caustic, here 
the amplitude for the single chord determined by 6_ and 6+ becomes singular but the 
integral is regularized by the determinant in ( |A.6| ). 

It is easier to deal first with the integral over the square of the Wigner function in 
( p.l4| ). Introducing in ( |1.4| ) the amplitude of a chord contribution to the semiclassical 
pure state Wigner function as B 



A{x) 



TT 



det{J_,/+}| 



(A.7) 
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and averaging over oscillations, i. e. taking cos^ 



i we immediately obtain 



trp" 



1. 



(A.8) 



(27r)' {27iy 

Given that the pure state condition is verified semiclassically, ( |A.8| ) indirectly should 
imply that, also tr p = 1. However, this is not the result of the direct calculation of 



tr p 



dxW(: 



1 



J d^_d^+ |det{I_,I+}|i cos 



H 



(A.9) 



(47rV27rfi)' 

This integral is dominated by the region where 6'+ is near 6-, even though the 
determinant in the r.h.s. of ([A .91 ) cancels for ^4. = 6-, because there the hessian 
determinant of S{9-,9^) also tends to zero. It is first necessary to estimate the area, 
5", for small chords. Let us return to the interpretation of a given angle, 9^, as the time 
along a trajectory on the torus, generated by choosing the respective action variable, I^, 
as the Hamiltonian. If we transfer the origin to the torus, at so that = = 0, 
the action for the chord to is just 

3(9^,9^)= f\{t)q{t) dt-ip+-q+. (A.IO) 

Hence, the second derivative of the chord action is simply 

2 



so that, as 6*4 



391 

9^ 



■1+ 

we obtain 



(A.ll) 



1 



X, A i_ 



(A.12) 



Therefore the hessian determinant for the stationary phase along 9^ 



9^ is just 



det 



Met{/_,/+}|. 



(A.13) 



As noted previously, the Poisson bracket matrix has / x / dimensions, but the power of 
its preceding factor in ( [A.ltj| ) is not the same as in ( |A.4| ), so that finally, 

trp^ . (A.14) 

The factor previously obtained using noncanonical coordinates on the torus, in 

for / 



1, was y2/3, instead of ([A.14| ). The important point is not the precise factor, 
but the fact that this depends neither on Planck's constant, nor on the classical system, 
except for the number of degrees of freedom. There is no point in trying to divide 
the semiclassical Wigner function by the normalization factor ( |A.14| ), because this will 
be inaccurate locally, inside the torus, as well as destroying the indirect normalization 
( |A.8| ). The failure of the normalization of the simple semiclassical Wigner function 
is a consequence of its (integrable) singularity precisely in the neighbourhood of the 
classical manifold, which is the region that dominates the normalization integral ( [A.9| ). 
In contrast, the integral for tr depends on the entire oscillatory domain of the Wigner 
function, so that indirect semiclassical normalization is exact. 



Decoherence of Wigner Functions 



23 



What is the effect of evolution on the normahzation? Evidently, both ( |A.8| ) and 
( |A.14| ) are stable with respect to unitary evolution. Because the normalization integral 
is dominated by the region of small chords, the decoherence distance functional ( |4.3| ) 
cancels for the stationary phase condition, 6^ = 6^, in ( [A .91 ), so that there is no effect of 
decoherence on the direct normalization. The evolving semiclassical Wigner function is 
given correctly by ( ^4.2| ) while the semiclassical approximation remains valid, i. e. until 
the amplitude becomes so concentrated on the classical manifold that the assumption 
of its smoothness is no loner valid. It follows that the remarkably simple expression 

f Af) AO 



in terms of the decoherence distance functional ( [4.3|) over the pair of orbits that end up 
on 6^ and is valid within the present semiclassical theory. 



